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This  study  conducts  experimental  investigation  and  numerical  analysis  for  one-stage  thermoelectric 
cooler  (TEC)  considering  Thomson  effect.  Three  Seebeck  coefficient  models  are  applied  to  numerically 
and  experimentally  study  the  Thomson  effect  on  TEC.  Results  show  that  higher  current,  higher  hot  side 
temperature,  or  lower  heat  load  can  increase  the  temperature  difference  between  the  cold  and  hot  sides. 
Opposite  trends  are  found  for  COP.  Specific  current  should  be  chosen  as  the  upper  threshold  in  thermo¬ 
electric  cooler  design.  The  cooling  performance  can  improve  when  the  Thomson  heat  maintains  positive. 
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1.  Introduction 

Development  of  high  density,  high  power,  and  high  efficiency  of 
electronic  components  enhances  the  system  efficiency  but  accom¬ 
panies  with  heat  dissipation  problems.  Dramatically  increasing  heat 
flux  may  result  in  the  failure  of  electronic  component  or  the  decrease 
in  efficiency  and  reliability  of  system.  Therefore,  an  efficient  cooling 
method  is  highly  desirable  to  precision  electronic  components. 

Thermoelectric  cooler  (TEC)  is  one  of  the  common  cooling 
methods  with  special  features  and  potential.  It  can  be  manufac¬ 
tured  from  millimeter  scale  to  micrometer  scale  and  even  can  be 
made  inside  the  electronic  component  to  do  thermal  management 
directly.  This  suits  well  for  the  limited  installation  space  inside 
electric  equipment  in  practical  application.  It’s  an  active  cooling 
method,  but  there  is  no  working  fluid  and  mechanical  operation  in¬ 
side.  Therefore,  its  silence  and  no  vibration  while  working  won’t 
cause  physical  effect  on  electronic  components.  Moreover,  due  to 
its  direct  conversion  between  electrical  energy  and  thermal  energy, 
thermal  management  is  more  efficient  than  other  electronic  cool¬ 
ing  devices  and  heat  can  be  controlled  precisely. 

Thermoelectric  effects  involve  the  conversion  between  electri¬ 
cal  energy  and  thermal  energy.  Three  main  findings  of  thermoelec¬ 
tric  effects,  Seebeck  effect,  Peltier  effect,  and  Thomson  effect,  can 
be  traced  back  to  the  early  1 9th  century.  Seebeck  effect  was  found 
that  in  a  closed  circuit  made  from  two  dissimilar  metals,  there 
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would  generate  an  electrical  current  with  the  junctions  at  different 
temperatures.  Seebeck  effect  can  occur  in  a  circuit  composed  of 
any  substance,  but  the  magnitude  of  the  Seebeck  coefficient  is 
dependent  on  the  material.  Peltier  effect  can  be  viewed  as  the  in¬ 
verse  effect  of  Seebeck  effect.  The  heat  absorbed  is  defined  as  the 
Peltier  heat.  Thomson  effect  was  observed  that  if  an  electrical  cur¬ 
rent  was  offered  to  a  single  conductor  with  a  temperature  differ¬ 
ence  between  two  endpoints,  material  will  emit  or  absorb  heat. 
The  combination  of  Seebeck  effect,  Peltier  effect  and  Thomson  ef¬ 
fect  is  called  “Kelvin  relationships”.  A  property,  “figure-of-merit 
(Z)”,  is  then  defined  by  scientists  to  represent  characteristics  of 
these  three  properties.  The  expression  of  figure-of-merit  is 


where  a  is  Seebeck  coefficient,  a  is  electrical  conductivity  and  k  is 
thermal  conductivity. 

Theories  of  thermoelectric  effects  were  fundamentally  devel¬ 
oped  by  scientists  [1-4].  Chen  et  al.  [5]  discussed  the  influence  of 
Thomson  effect  on  the  thermoelectric  generator  and  provided  the 
new  relationships  considering  Thomson  effect  to  express  the  max¬ 
imum  power  output  and  the  maximum  efficiency.  Yamanashi  [6] 
proposed  to  optimize  the  design  of  the  thermoelectric  cooler  by 
the  dimensionless  entropy  flow  equations.  Chua  et  al.  [7]  studied 
the  relationship  between  temperature  and  entropy.  The  entropy 
density  was  applied  to  explain  the  capacity  of  thermoelectric  cool¬ 
ing  and  thermoelectric  generation.  Xuan  et  al.  [8-13]  adopted  fun¬ 
damental  thermodynamics  and  temperature-entropy  diagrams  to 
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Nomenclature 

4 

cross  section  area  of  one  semiconductor  pellet 

AT 

temperature  difference  between  hot  side  and  cold  side 

COP 

coefficient  of  performance 

CSM 

constant  Seebeck  coefficient  model 

Greek  symbols 

i 

variable 

a 

Seebeck  coefficient 

I 

current 

K 

heat  conductivity 

J 

current  density 

Pt 

Thomson  heat 

LLS 

log-linear  Seebeck  coefficient  model 

P 

electrical  resistivity 

LSM 

linear  Seebeck  coefficient  model 

G 

electrical  conductivity 

lx 

length  of  one  semiconductor  pellet 

T 

Thomson  coefficient 

N 

number  of  P-N  pair 

Ptec 

thermoelectric  cooler  power 

Subscripts 

Qc 

cold  side  heat  load 

c 

cold  side 

Qu 

hot  side  heat  load 

max 

maximum 

R2 

coefficient  of  determination 

opt 

optimum 

T 

temperature 

rev 

reverse 

Tc 

cold  side  temperature 

tau 

positive  Thomson  heat 

Tu 

hot  side  temperature 

TEC 

thermoelectric  cooler 

VjEC 

thermoelectric  cooler  voltage 

u 

hot  side 

X 

the  x-axis  in  Cartesian  coordinate  system 

X 

x-direction 

z 

Figure  of  merit 

analyze  the  internal  heat  distribution  of  thermoelectric  cooler  at 
steady  state  and  when  the  interface  layer  is  considered.  Multi¬ 
stage  thermoelectric  coolers  were  further  investigated  and  the  po¬ 
lar  characteristics  of  the  coefficient  of  performance  were  used  to 
optimize  the  design.  Yang  and  Chen  [14]  analyzed  the  cooling 
capacity  of  one-  and  two-stage  thermoelectric  microcoolers  with¬ 
out  considering  Thomson  effect  and  pointed  that  the  techniques 
of  integrated  circuit  and  microelectromechanical  manufacturing 
are  two  primary  skills  for  thermoelectric  cooling  development. 
Chen  et  al.  [15]  applied  thermoelectric  cooler  to  high  power 
light-emitting  diode  (LED)  by  microfabrication  technique  and 
flip-chip  assembly  process,  and  analyzed  the  effect  of  varied  semi¬ 
conductor  pairs  on  decreasing  the  internal  heat  resistivity  of  LED. 
Hodes  [16]  simulated  one-dimensional  temperature  distribution 
of  thermoelectric  modules.  Different  from  most  previous  investiga¬ 
tions,  the  controlled  side  and  uncontrolled  side  were  defined  as 
boundary  conditions.  The  controlled  side  had  three  operating 
forms  with  varied  current  supplied:  cooling,  heating  and  neutral 
forms.  Huang  et  al.  [17]  studied  the  influence  of  Thomson  effect 
by  simulating  the  temperature  distribution  and  coefficient  of  per¬ 
formance  of  one-stage  thermoelectric  cooler  with  known  cold  side 
and  hot  side  temperatures.  Results  showed  that  Thomson  heat 
could  lower  the  effect  of  Joule  heat  at  the  cold  side.  The  methods 
to  improve  the  coefficient  of  performance  should  consider  not  only 
a  better  figure-of-merit  but  also  the  influence  of  Thomson  effect. 
Kraftmakher  [18]  proposed  how  to  measure  the  Seebeck  coefficient 
and  the  Peltier  coefficient  of  thermoelectric  cooler  with  simple 
experimental  methods.  Lee  and  Kim  [19]  conducted  three-dimen¬ 
sional  numerical  analysis  o  predict  the  optimal  performance  for 
thermoelectric  micro-cooler.  Results  showed  that  the  performance 
of  thermoelectric  cooler  can  be  improved  more  for  thick  thermo¬ 
electric  elements  at  small  temperature  difference  and  small  cur¬ 
rent.  The  key  factors  to  evaluate  the  performance  of  the 
thermoelectric  microcooler  were  more  pairs  of  semiconductor  or 
smaller  cross-sectional  area  of  semiconductor.  Chang  et  al.  [20] 
established  a  theoretical  model  of  the  thermal  analogy  network 
for  thermoelectric  air-cooling  module.  The  influence  of  thermal 
resistances  was  analyzed  and  the  appropriate  operating  range 
was  predicted.  Yamashita  [21]  investigated  the  resultant  Seebeck 
coefficient  which  is  derived  from  the  temperature  dependence  of 
the  intrinsic  Seebeck  coefficient  by  considering  the  Thomson  effect. 
It  is  found  that  both  the  resultant  Seebeck  coefficient  and  the 


resultant  figure-of-merit  were  larger  than  the  intrinsic  properties 
and  the  figure-of-merit  was  affected  dramatically  at  the  interfaces 
when  the  intrinsic  Seebeck  coefficient  varied  nonlinearly  with 
temperature.  Pichanusakorn  and  Bandaru  [22]  maximized  the 
power  factor  (a2 a)  to  find  the  optimal  Seebeck  coefficient.  Results 
showed  that  for  any  material,  the  maximum  power  factor  at  a  gi¬ 
ven  temperature  is  expected  when  the  Seebeck  coefficient  in  the 
range  of  130-187  pV/K. 

In  previous  investigations  of  thermoelectric  cooler,  constant  hot 
and  cold  side  temperatures  of  thermoelectric  module  were  often 
chosen  as  boundary  conditions.  However,  like  Hodes  [16],  the  hot 
side  temperature  and  the  cold  side  heat  load  are  controlled  in  this 
study.  The  cold  side  heat  load  is  controlled  in  order  to  meet  the 
known  allowable  power  input  for  the  electronic  components  in 
practical  application.  The  hot  side  temperature  of  thermoelectric 
module  remains  fixed  to  easily  observe  the  cold  side  temperature 
variation.  Therefore,  with  a  variety  of  input  current,  the  heat  load 
of  the  cold  side,  and  the  hot  side  temperature  controlled  by  the 
blower  system,  the  cold  side  temperature  and  the  coefficient  of 
performance  of  one-stage  thermoelectric  cooler  are  obtained 
through  the  experiment  and  numerical  simulation. 

On  the  other  hand,  not  only  Joule  heat  but  also  Thomson  heat  is 
taken  into  account  in  the  internal  thermal  analysis  of  thermoelec¬ 
tric  cooler.  In  order  to  discuss  the  effect  of  Thomson  heat  on  tem¬ 
perature  prediction  and  the  interior  heat  transfer  mechanism, 
three  different  Seebeck  coefficient  models,  constant  Seebeck  model 
(CSM),  second-order  linear  Seebeck  model  (LSM),  and  first-order 
log-linear  Seebeck  model  (LLS),  are  examined  through  experimen¬ 
tal  investigation  and  numerical  simulation.  CSM  indicates  a  zero 
Thomson  coefficient.  LSM  means  the  Thomson  coefficient  is  a  func¬ 
tion  of  temperature,  LLS  represents  a  constant  Thomson  coefficient. 
The  best  Seebeck  coefficient  model  is  determined  by  comparing 
simulation  results  with  experimental  results.  Thermal  analysis 
considering  Thomson  heat  is  then  simulated  with  the  best  Seebeck 
model,  and  the  influence  of  Thomson  heat,  Fourier  heat  and  Joule 
heat  on  one-stage  thermoelectric  cooler  is  analyzed. 


2.  General  model  of  one-stage  thermoelectric  cooler 

A  thermoelectric  cooler  is  composed  of  several  semiconductor 
pairs  which  are  sandwiched  between  two  ceramic  layers,  as  shown 
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Uncontrolled  Side  Interface,  usi 
Controlled  Side  Interface,  csi 


Fig.  1.  Schematic  diagram  of  one  pair  of  semiconductors. 


in  Fig.  1.  A  pair  of  semiconductors  in  thermoelectric  cooler  is  cho¬ 
sen  as  the  micro-system  for  analysis.  The  x-axis  of  the  thermoelec¬ 
tric  cooler  is  defined  positive  from  down  side  to  upper  side,  i.e.  the 
flowing  direction  of  current  in  p-type  is  positive.  The  down  side 
where  the  heat  absorption  occurs  is  controlled  side,  the  cold  side 
of  thermoelectric  cooler.  The  upper  side  where  the  heat  release  oc¬ 
curs  is  uncontrolled  side,  the  hot  side  of  thermoelectric  cooler. 
Properties  of  p-type  and  n-type  are  assumed  to  be  equal  except 
the  Seebeck  coefficient  and  electric  current  are  negative  in  n-type. 
Thus  the  temperature  distributions  in  p-type  and  in  n-type  are  the 
same.  Governing  equation  of  one-stage  thermoelectric  cooler  is  de¬ 
rived  according  to  irreversible  thermodynamics,  thermoelectric 
relations  and  Onsager’s  Reciprocal  Relations  [1-4].  One  can  obtain 
the  three-dimensional  governing  equation 

kV2T  —  xj  ■  VI  +  pj2  =  0,  (2) 

where  k  is  thermal  conductivity,  t  is  Thomson  coefficient,  J  is  de¬ 
fined  as  the  electric  current  density,  and  electrical  resistivity  p  is 
the  reciprocal  of  electrical  conductivity  (p  =  l/cr).  The  first  term  is 
the  heat  resulted  from  Fourier  conduction,  the  second  term  and 
the  last  term  are  attributed  to  heat  generation  from  the  Thomson 
heat  (px)  and  the  Joule  heat.Different  from  previous  investigations 
using  constant  hot  and  cold  side  temperatures  as  boundary  condi¬ 
tions,  the  hot  side  temperature  ( Tu )  and  the  cold  side  heat  load 
(Qc)  are  controlled  in  this  study.  The  boundary  equations  are 

TU  =  Tu,  (3) 


QU  =  Qc  (4) 

The  heat  balance  at  the  hot  side  boundary  is 

Qu  =  2MaU/u-2NKAp^Ux,  (5) 

where  N  is  the  number  of  pair  of  semiconductors;  Ap  is  the  cross- 
sectional  area  of  p-type  semiconductor.  The  first  term  indicates 
the  heat  resulted  from  Peltier  cooling  effect  and  the  second  term 
is  from  the  Fourier  conduction.  The  total  power  of  one-stage  ther¬ 
moelectric  cooler  can  be  obtained  from  the  heat  loads, 


Ptec  =  Qu  -  Qc  =  2M«UxTu  -  2 NkA^ 


-Qc- 


x=Lx 


(6) 


Therefore,  the  input  voltage  of  one-stage  thermoelectric  cooler  is 


Vtec  =  2Noc\x=lTu 


2N-— 

Jdx 


Qc 
I  ' 


(7) 


Finally,  the  coefficient  of  performance  (COP)  of  test  module  is  de¬ 
fined  as  the  cooling  rate  per  unit  power  input  of  test  module, 


COP  =  = _ — _ 

Ptec  2MaU/U  “  2N^r*U  “  Q^ ' 

Only  when  AT  >  0,  COP  has  significance. 


(8) 


3.  Numerical  methods 


In  this  study  Thomson  heat  (pT)  is  taken  into  account  in  order  to 
discuss  its  effect  on  temperature  prediction  and  the  internal  heat 
transfer  mechanism  of  one-stage  thermoelectric  cooler.  Thomson 
heat  and  Thomson  coefficient  (t)  can  be  determined  by  the  See¬ 
beck  coefficient.  Therefore,  three  different  Seebeck  coefficient 
models,  constant  Seebeck  model  (CSM),  second-order  linear  See¬ 
beck  model  (LSM),  and  first-order  log-linear  Seebeck  model  (LLS), 
are  examined  for  suitability  and  accuracy.  The  mathematical  forms 
of  these  three  models  are  listed  in  Table  1.  In  CSM  model,  Seebeck 
coefficient  is  a  constant  property,  which  is  mostly  used  in  thermo¬ 
electric  cooler  analysis.  The  constant  Seebeck  coefficient  leads  to  a 
zero  Thomson  Coefficient,  i.e.  the  Thomson  effect  is  neglected.  LSM 
model  and  LLS  model  are  both  temperature-dependent.  LSM  model 
indicates  that  both  Seebeck  coefficient  and  Thomson  coefficient  are 
the  function  of  temperature.  LLS  model  has  a  constant  Thomson 
coefficient.  The  coefficients  C0,  Ci  and  C2  in  these  three  models 
are  determined  through  experiment. 

Several  assumptions  are  made  to  simplify  the  simulation:  (1) 
The  test  module  operates  in  steady  state,  and  all  properties  are 
independent  of  time;  (2)  Electric  current  flows  along  single  direc¬ 
tion  without  scattering;  (3)  Energy  only  transmits  across  the  cold 
side  and  the  hot  side  of  the  test  module.  The  other  parts  are  as¬ 
sumed  to  be  in  adiabatic  condition,  i.e.  no  energy  transfer;  (4)  Heat 
only  flows  along  the  x-direction;  (5)  Heat  load  supplied  by  heating 
unit  is  completely  absorbed  by  the  test  module;  (6)  Thermal  con¬ 
ductivity  and  electrical  resistivity  are  viewed  as  constants  at  any 
temperature,  and  their  values  are  obtained  from  the  average  value 
of  measurement;  (8)  Seebeck  coefficient  at  different  temperatures 
is  measured  in  experiment.  The  data  of  Seebeck  coefficient  will  be 
applied  to  determine  CSM  model,  LSM  model  and  LLS  model;  (9) 
The  properties  of  p-type  are  equal  to  the  properties  of  n-type.  Only 
the  electric  current  and  the  Seebeck  coefficient  are  negative  in  n- 
type. 

3.1.  Second-order  linear  Seebeck  model  (LSM) 


Second-order  LSM  model  represents  that  the  Thomson  coeffi¬ 
cient  is  a  second-order  polynomial.  Hence  the  governing  equation 
becomes  a  non-linear  ordinary  differential  equation 


d2r 

dx2 


-  (2  c2r2  +  c,  T)  ~ +L?  =  o. 

K v  J  dx  K 


(9) 


The  governing  equation  is  then  resolved  numerically 


^2  j(t+l) 

dx2 


<t+i) 


J-(2C2T2  +  QTf)dT- 
kk  j  dx 


+&  =  o, 


(10) 


where  the  superscript  (t)  indicates  the  old  value;  (t+  1)  means  the 
new  value.Eq.  (10)  is  rewritten  with  finite  difference  method, 


Table  1 

Mathematical  form  of  Seebeck  coefficient  models  in  this 
study. 


Seebeck  models 

Formula 

2nd-order  LSM 

a  =  C0  +  C\T  +  C2T2 

CSM 

a  =  C0 

1  st-order  LLS 

a  =  C0  +  Ci  In  T 
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Ti+1  -  2Ti  +  Ti_,\ (t+1) 


J 

K 

=  0. 


Ax2 

[2C2(T(t))2  +  C,T(t) 


T  T  \  (f+L  j2 

1  ;+ 1  —  1  i  ^  x 


2Ax 


+' 


Tp 


(11) 


Eq.  ( 1 1 )  is  rearranged, 

_fpAx^  =  _  J-Ax[2C2  (Tm)2  +  C1T(t)]}r^,)  -  2Tit+V 

+  {1  +^Ax[2  C2(T(t))2  +  C1T(t)]}T(iV)- 

The  heat  balance  at  the  cold  side  boundary  is 

Qc  =  2NIajx=0T\x=0  -  2NkAp^ 


Q c  rr  ,  r  i  r  rrm2iT  K  f-3Ti-\-4Ti+i  -Ti+2 
2M  =  [C°  +  Cl  1(0)  +  C2I(0)  ]T,  -  J  - 2Ax - 

Therefore,  Eq.  (14)  is  rewritten  in  the  form  of  point  solution, 

Qc  I  _K_  MtQ)  _  t(Q  \ 
j(t+ 1)  _  2M  ~r~  2AxJ  i+1  1  i+2  / 

'■  _C0  +  C17f  +  C2(7f)2+2|;' 


(12) 

(13) 

(14) 

(15) 


3.2.  Constant  Seebeck  model  (CSM) 


Thomson  coefficient  (t)  is  zero  in  CSM.  The  governing  equation 
is 


^+^  =  0. 
dx2  k 


(16) 


The  equation  turns  into  a  second-order  linear  ordinary  differential 
equation  and  can  be  solved  analytically. 


3.3.  First-order  log-linear  Seebeck  model  (LSM) 


The  characteristic  in  lst-order  LSM  is  that  Thomson  coefficient 
is  a  constant  in  this  model.  Then  the  governing  equation  is  rewrit¬ 
ten  as 


d2r 

dx2 


c,“!r+— ■»- 

K  dx  K 


(17) 


The  governing  equation  also  becomes  a  linear  problem.  Therefore,  it 
can  be  either  solved  analytically  like  CSM  model  or  numerically  like 
LSM  model. 


pairs  of  semiconductors  and  two  interface  layers.  The  semiconduc¬ 
tor  material  is  Bismuth  Telluride  (BiTe).  Height  is  about  2  mm  and 
cross-sectional  area  is  1.4  mm2.  The  two  0.9  mm2  Alumina  ceramic 
interface  layers  (A1203)  are  thermally  conductive  but  electrically 
insulated.  Bismuth  Tin  (BiSn)  is  used  to  connect  the  interface  layer 
and  semiconductors.  In  between  semiconductors,  there  is  silicone, 
RTV  filled  to  prevent  from  vapor.  The  heat  spreader  is  used  to  dis¬ 
tribute  the  heat  uniformly.  It  is  made  of  brass  which  has  high  ther¬ 
mal  conductivity  and  is  tractable  and  dysoxidizable.  Three  grooves 
are  machined  on  the  spreader  to  set  the  thermocouples  without 
destroying  the  surface  of  thermoelectric  coolers.  The  insulating 
blanket  (cerablanket)  made  by  Thermal  Ceramic  Inc.  is  made  of 
ceramic  fiber  materials,  like  silica  (Si02)  and  Aluminum  oxide 
(A1203).  It  has  very  low  thermal  conductivity  about  0.07  W/m  K  be¬ 
low  200  °C  to  assure  the  adiabatic  condition  around  the  test  mod¬ 
ule.  The  fins  are  used  as  the  heat  sink  and  assembled  at  the  top  of 
test  module.  Totally  144  fins  with  1.45  mm2  rectangle  cross-sec¬ 
tional  area  and  45.5  mm  height  are  made  of  A1  6063  in  one  work 
piece  and  are  uniformly  set  on  the  base  plate  with  45  mm2  area 
and  4.5  mm  height.  A  thin  silicon  rubber  heater  (40  mm2)  is  used 
as  the  heating  element  to  meet  the  known  heat  load  of  electronic 
component  in  practice.  Since  the  electrical  resistance  of  the  heater 
is  fixed,  the  heat  load  is  then  controlled  by  the  input  current.  In 
addition,  a  complementary  heater  is  used  to  make  sure  that  the 
heat  will  not  dissipate  from  the  bottom.  All  junctions  between 
components  are  filled  with  thermal  compound,  Arctic  Silver  5  type, 
made  by  Arctic  Silver  Incorporated  Company.  The  average  particle 
size  of  the  thermal  compound  is  small  enough  to  fill  the  cavities  at 
the  contact  surfaces  in  order  to  reduce  the  contact  thermal 
resistance. 

The  data  acquisition  system  includes  thermocouple,  data  acqui¬ 
sition  equipment  (DAQJ  and  personal  computer.  The  type  T  ther¬ 
mocouple  is  used  to  measure  the  hot  side  temperature  and  cold 
side  temperature  of  the  test  module.  The  data  acquisition  equip¬ 
ment,  NetDAQ  2640A,  made  by  Fluke  Corporation  Company  in¬ 
cludes  a  universal  input  module  with  20  analog  channels,  10 
computed  channels  and  8  alarm  channels.  It  can  measure  DC  and 
AC  voltage,  resistance  and  temperature. 

As  shown  in  Fig.  4,  surroundings  control  system  composed  of 
fan,  wind  tunnel,  thermostat,  and  settling  chamber  is  utilized  to 
maintain  the  desired  hot  side  temperature  of  test  module.  The 
fan  drives  the  air  to  flow  along  the  wind  tunnel  and  passes  through 
the  fin  at  the  top  of  test  section  to  take  the  heat  away  from  the  test 
section.  The  settling  chamber  made  of  a  porous  aluminum  board  is 
installed  at  the  entrance  of  the  wind  tunnel  in  order  to  assure  that 
the  air  can  flow  more  stably.  The  thermostat  controls  the  air 
properties. 


4.  Experimental  methods 

4.1.  Test  facility 

As  shown  in  Fig.  2,  the  main  components  of  the  test  apparatus 
consist  of  test  section,  power  supply,  data  acquisition  system, 
and  surroundings  control  system.  The  linear  DC  power  supply 
made  by  the  Good  Will  Instrument  Co.  Ltd.  provides  the  required 
voltage  and  current  for  test  section  and  the  forced  convection  sec¬ 
tion  in  surroundings  control  system. 

Fig.  3  shows  the  detailed  construction  of  the  test  section.  It 
mainly  includes  test  module,  heat  sink  unit,  cerablanket,  and  heat¬ 
ing  unit.  The  test  module  is  composed  of  one  thermoelectric  cooler 
and  two  heat  spreaders.  Thermoelectric  cooler  is  TEC1 -127.08, 
made  by  Tande  Energy  and  Temperature  Associates  Pty.  Co.  Thick¬ 
ness  is  3.8  mm,  working  area  is  40  mm2,  and  the  operation  current 
range  is  0-8  A.  The  thermoelectric  cooler  is  constructed  from  125 


4.2.  Experimental  procedure 

Experiments  are  divided  into  pre-experiment  and  main  experi¬ 
ment.  Three  important  properties  of  thermoelectric  cooler,  thermal 
conductivity,  electrical  resistivity  and  Seebeck  coefficient,  are  mea¬ 
sured  in  pre-experiment.  The  hot  side  temperature,  cold  side  tem¬ 
perature,  voltage  and  current  of  one-stage  thermoelectric  cooler 
are  measured  in  main  experiment. 

In  pre-experiment,  there  is  no  need  to  supply  current  to  the 
thermoelectric  cooler.  The  temperature  difference  of  the  thermo¬ 
electric  cooler  is  measured  with  a  given  heat  load.  One-dimen¬ 
sional  heat  conduction  is  assumed  and  the  thermal  conductivity 
at  different  temperatures  can  be  calculated  by  Fourier’s  law.  The 
electrical  resistance  is  measured  at  different  temperatures.  By  con¬ 
trolling  the  temperature  difference  at  2.5  K  ±  0.3  K,  the  voltage 
generating  by  the  temperature  difference  is  measured  at  different 
temperatures.  The  Seebeck  coefficient  is  then  determined  by  the 
ratio  of  the  voltage  and  the  temperature  difference. 
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Fig.  2.  Experimental  facility. 


Fig.  3.  Construction  of  test  section. 


Fig.  4.  Construction  of  surroundings  control  system. 


In  main  experiment,  three  hot  side  temperatures,  323.2  K, 
333.2  K,  and  343.2  K  are  selected.  Three  input  heat  loads,  5.5  W, 


12.3  W  and  22.0  W  are  determined  by  the  input  current  of  heating 
element  with  0.5  A,  0.75  A  and  1.0  A.  Seven  input  currents  of  test 
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module,  1.0  A,  1.5  A,  2.0  A,  2.5  A,  3.0  A,  3.5  A  and  4.0  A  are  exam¬ 
ined.  The  cold  side  heat  load  is  controlled  by  the  heating  element 
whose  current  is  supplied  by  the  power  supply.  The  hot  side  tem¬ 
perature  is  controlled  by  the  forced  convection  section.  Once  the 
change  of  the  hot  side  temperature  is  in  ±0.6  K  during  900  s,  it 
reaches  steady  state.  The  cold  side  temperature  and  the  voltage 
of  test  module  are  then  collected  and  stored. 

5.  Results  and  discussion 

5.1.  Properties  of  thermoelectric  cooler 

Thermal  conductivity  and  electrical  resistivity  of  single  thermo¬ 
electric  cooler  are  first  measured  at  different  temperatures.  Results 
show  that  thermal  conductivity  decreases  but  electrical  resistivity 
increases  with  increasing  temperature.  However,  both  changes  are 
small  enough  to  be  neglected.  Thus  thermal  conductivity  and  elec¬ 
trical  resistivity  can  be  assumed  constant.  The  properties  of  ther¬ 
moelectric  cooler  are  listed  in  Table  2. 

Fig.  5  shows  Seebeck  coefficient  is  a  function  of  temperature 
and  increases  with  increasing  temperature.  Three  Seebeck  coeffi¬ 
cient  models  (CSM,  LSM  and  LLS)  are  examined  for  suitability 
and  accuracy.  The  average  value  of  all  measured  data  is  used  as 
the  constant  Seebeck  coefficient  in  CSM  model.  The  formulas  of 
LSM  model  and  LLS  model  are  obtained  by  regression.  On  the  other 
hand,  Thomson  coefficient  is  defined  as 


Fig.  6  shows  the  Thomson  coefficient  for  three  models.  CSM  model 
leads  to  a  zero,  LSM  model  is  temperature-dependent,  and  LLS  mod¬ 
el  has  a  constant  value. 

5.2.  Experimental  results 

For  one-stage  thermoelectric  cooler  with  fixed  hot  side  temper¬ 
ature  (Tu)  and  the  fixed  cold  side  heat  load  (QJ,  the  temperature 
difference  (AT)  and  COP  are  obtained  by  changing  the  input  cur¬ 
rent  of  test  module  (/).  Fig.  7  shows  the  distributions  of  tempera¬ 
ture  difference  with  increasing  current  for  different  Tu  and 
different  Qc.  Due  to  the  insufficiency  of  the  surroundings  control 
system  when  a  low  current  and  a  high  cold  side  heat  load  are  ap¬ 
plied,  the  desired  hot  side  temperature  cannot  be  attained  so  that 
there  is  no  measured  data.  Results  show  that  temperature  differ¬ 
ence  increases  with  increasing  current.  Since  the  hot  side  temper¬ 
ature  is  fixed,  this  means  the  cold  side  temperature  goes  down 
with  increasing  current.  At  the  same  current,  temperature  differ¬ 
ence  decreases  with  increasing  heat  load.  It  is  observed  that  the 
slope  of  temperature  difference  decreases  with  increasing  current 
and  is  expected  to  approximate  to  zero  at  a  specific  current  which 
means  there  will  exist  a  maximum  temperature  difference.  Regard¬ 
ing  the  effect  of  hot  side  temperature,  overall  the  trends  are  found 
similar.  To  emphasize  the  difference  in  magnitude,  Fig.  8  shows  the 
distributions  of  temperature  difference  with  increasing  current  for 
different  Tu  at  Qc  =  12.3  W.  At  the  same  current,  temperature  differ- 


Table  2 

Properties  of  one-stage  TEC  in  this  study. 


Properties 

Values 

Number  of  P-N  pair  (N) 

Length  of  one  semiconductor  pellet  (Lx) 

Cross-sectional  area  of  one  semiconductor  pellet  (Ap) 
Thermal  conductivity  of  one  semiconductor  pellet  (k) 
Electrical  resistance  of  one  semiconductor  pellet  (R) 
Electrical  resistivity  of  one  semiconductor  pellet  (p) 

125 

2.0  x  10“3m 

1.96  x  10“6  m2 

2.571  W  m-1  K1 

1.80  Q 

7.056  x  10-6  Q  m 

Fig.  6.  Comparison  of  Thomson  coefficient  with  different  Seebeck  coefficient 
models. 

ence  slightly  increases  with  increasing  hot  side  temperature.  How¬ 
ever,  comparing  with  the  effect  of  heat  load,  the  change  of 
temperature  difference  resulted  from  the  different  hot  side  tem¬ 
perature  is  minor. 
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Fig.  7.  Temperature  difference  at  different  hot  side  temperatures  with  different 
cold  side  heat  loads. 


Fig.  8.  The  variations  of  temperature  difference  at  different  hot  side  temperatures 
with  the  cold  side  heat  load  12.3  W. 


Fig.  9  shows  the  distributions  of  COP  with  increasing  current  for 
different  Tu  with  different  Qc-  The  COP  has  no  meaning  when  the 
test  module  loses  its  cooling  capacity  (AT  <  0).  Therefore,  the 
meaningless  data  will  not  be  shown  in  the  figure.  The  trend  of 
COP  is  opposites  to  that  of  AT,  the  COP  decreases  with  increasing 
current.  According  to  the  definition  of  COP,  Eq.  (8),  the  COP  de¬ 
creases  when  increasing  the  current  of  test  module  with  a  fixed 
cold  side  heat  load.  Moreover,  the  COP  decreases  quickly  at  higher 
current  and  almost  approximates  to  zero.  The  maximum  COP  ap¬ 
pears  when  the  TEC  current  is  small  enough  to  make  the  tempera¬ 
ture  differences  approximate  to  zero.  With  fixed  current,  adding 
heat  load  makes  COP  larger  and  temperature  difference  smaller. 
Fig.  10  shows  the  difference  of  hot  side  temperature  has  relatively 
small  effect  on  COP  and  can  almost  be  neglected. 

5.3.  Numerical  simulation 

Higher  current  values  are  examined  in  order  to  find  the  maxi¬ 
mum  temperature  difference  in  numerical  simulation.  Thomson 
heat  is  taken  into  account  in  order  to  discuss  its  effect  on  the  inter¬ 
nal  temperature  prediction  and  heat  transfer  mechanism  for  one- 
stage  thermoelectric  cooler. 

In  Fig.  11,  the  temperature  differences  for  three  different  See- 
beck  models  are  simulated  at  the  hot  side  temperature  343.2  K 
with  the  cold  side  heat  load  12.3  W.  Results  show  that  there  does 
exist  a  maximum  temperature  difference  (A Tmax)  for  each  model 
at  a  specific  current  which  is  defined  as  the  optimum  current,  Iopt. 
Larger  than  this  current,  the  heat  generation  will  generate  faster 
than  the  Peltier  cooling  heat,  so  the  temperature  difference  starts 


4882 


C.-Y.  Du,  C.-D.  Wen  /  International  Journal  of  Heat  and  Mass  Transfer  54  (2011)  4875-4884 


Fig.  10.  The  variations  of  COP  at  different  hot  side  temperatures  with  the  cold  side 
heat  load  12.3  W. 


Fig.  11.  Comparison  of  temperature  difference  simulated  by  CSM  model,  LSM 
model  and  LLS  model. 


to  decrease.  Therefore,  in  thermoelectric  cooler  design  the  input 
current  must  be  smaller  than  the  optimal  current  in  order  to 


achieve  cooling.  In  addition,  the  effect  of  Seebeck  model  becomes 
apparent  at  higher  current  range.  It  is  necessary  to  examine  their 
suitability  and  accuracy  while  being  applied  in  simulation. 

Figs.  12  and  13  show  the  comparison  of  numerical  results  and 
experimental  data  for  temperature  difference  and  COP  at  the  hot 
side  temperature  343.2  I<  with  the  cold  side  heat  load  12.3  W.  From 
the  figures,  LSM  model  shows  the  best  agreement  with  the  exper¬ 
imental  data,  CSM  model  in  second,  and  LLS  model  in  last.  In  LSM 
model,  the  Seebeck  coefficient  and  the  Thomson  coefficient  are 
both  dependent  on  temperature.  The  poor  results  of  LLS  model 
may  be  due  to  the  over-predicting  by  constant  Thomson  coeffi¬ 
cient.  Overall,  LSM  model  is  found  to  be  the  best  Seebeck  coeffi¬ 
cient  model  in  this  study. 

The  internal  temperature  prediction  and  heat  transfer  mecha¬ 
nism  for  one-stage  thermoelectric  cooler  is  simulated  with  the  best 
Seebeck  model,  LSM  model.  At  a  given  hot  side  temperature 
333.2  K  with  the  cold  side  heat  load  12.3  W,  the  internal  tempera¬ 
ture  distributions  with  current  change  (0-12  A)  is  shown  in  Fig.  14. 
When  the  current  is  zero,  only  Fourier  conduction  occurs  so  that 
the  temperature  at  the  cold  side  near  the  heat  source  is  higher  than 
the  hot  side  temperature  near  surroundings.  Thermoelectric  effects 
occur  after  providing  the  current  into  the  test  module  and  cause 
the  drop  of  cold  side  temperature.  The  internal  temperature  will 
be  larger  than  the  hot  side  temperature  beyond  a  specific  current 
which  is  defined  as  reverse  current,  Jrev,  where  the  positive  temper¬ 
ature  gradient  at  the  hot  side  starts  to  turn  into  negative.  The  heat 
generation  inside  will  increase  faster  than  the  cooling  from  the  hot 
side.  Therefore,  for  practical  use,  input  current  should  be  lower 
than  lrev  to  prevent  from  high  internal  temperature  which  may  re¬ 
sult  in  the  failure  of  test  module.  As  defined  earlier,  the  optimum 
current,  Iopt ,  will  cause  the  maximum  temperature  difference,  i.e. 
the  cold  side  temperature  is  the  lowest.  According  to  the  results, 


Fig.  12.  Comparison  of  temperature  difference  simulated  by  three  Seebeck  coef¬ 
ficient  models  at  the  hot  side  temperature  343.2  K  with  the  cold  side  heat  load 
12.3  W. 
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Current  of  TEM  (A) 

Fig.  13.  Comparison  of  COP  simulated  by  three  Seebeck  coefficient  models  at  the 
hot  side  temperature  343.2  K  with  the  cold  side  heat  load  12.3  W. 


0.0  0.2  0.4  0.6  0.8  1.0 


x/Lt 

Fig.  14.  The  internal  temperature  distribution  versus  different  currents  at  the  hot 
side  temperature  333.2  K  with  the  cold  side  heat  load  12.3  W  using  LSM  model. 

/opt  =8.92  A  and  Jrev  =  3.66  A,  lrev  should  be  chosen  as  the  upper 
threshold  in  thermoelectric  cooler  design. 


x/L, 

Fig.  15.  Comparison  of  internal  temperature  distribution  between  LSM  model  and 
CSM  model. 


The  internal  temperature  distribution  of  LSM  model  considering 
Thomson  heat  and  CSM  model  neglecting  Thomson  heat  is  shown 


Fig.  16.  Distribution  of  Thomson  heat  at  the  hot  side  temperature  333.2  K  and  with 
the  cold  side  heat  load  12.3  W  using  LSM  model. 
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in  Fig.  15.  As  discussed  earlier,  the  cold  side  temperature  predic¬ 
tion  in  LSM  model  is  better  than  in  CSM  model  by  comparing  with 
experimental  data.  Results  here  show  the  temperature  simulated 
in  LSM  model  is  obviously  lower  than  in  CSM  model.  This  indicates 
that  applying  the  CSM  model  is  a  more  conservative  method  to 
estimate  the  maximum  temperature  difference  of  TEC,  while  the 
actual  temperature  difference  is  larger  in  application. 

The  internal  thermal  mechanism  includes  Fourier  heat  conduc¬ 
tion,  Joule  heat  and  Thomson  heat.  The  sum  of  Thomson  heat  and 
Joule  heat  plays  the  role  of  heat  generation  inside  TEC.  The  distri¬ 
bution  of  Thomson  heat  is  shown  in  Fig.  16.  When  the  Thomson 
heat  is  positive,  it  plays  the  role  of  cooling  in  the  interior  to  reduce 
the  effect  of  Joule  heat  and  Fourier  conduction  and  enhance  the 
cooling  capacity  of  the  TEC.  If  the  Thomson  heat  at  any  position  in¬ 
side  turns  into  negative,  where  the  current  is  defined  as  Itaw  the 
Thomson  heat  becomes  part  of  the  heat  generation  and  causes 
the  internal  temperature  to  rise  quickly.  Therefore,  the  cooling  per¬ 
formance  can  improve  when  the  Thomson  heat  maintains  positive. 
Comparing  the  value  of  Joule  heat  and  Thomson  heat,  Joule  heat  is 
larger  than  Thomson  heat.  Hence  the  Joule  heat  dominates  in  the 
internal  thermal  mechanism,  but  there  is  a  need  to  achieve  more 
precise  temperature  control  by  considering  Thomson  heat. 

6.  Conclusions 

The  purpose  of  this  research  is  to  analyze  the  performance  of 
one-stage  thermoelectric  coolers  by  experiment  and  numerical 
simulation.  The  influence  of  Thomson  effect  is  also  discussed. 
Key  findings  from  the  study  are  summarized  below. 

1.  Temperature  difference  increases  with  increasing  current  and 
decreases  with  increasing  heat  load.  Temperature  difference 
slightly  increases  with  increasing  hot  side  temperature,  but 
the  effect  of  hot  side  temperature  is  minor. 

2.  COP  decreases  with  increasing  current  and  approximates  to 
zero  at  a  high  current.  The  maximum  COP  appears  when  the 
TEC  current  is  small  enough  to  make  the  temperature  differ¬ 
ences  approximate  to  zero.  Adding  heat  load  makes  COP  larger 
and  temperature  difference  smaller.  The  hot  side  temperature 
has  relatively  small  effect  on  COP  and  can  almost  be  neglected. 

3.  There  exists  a  maximum  temperature  difference  at  high  cur¬ 
rent.  The  effect  of  Seebeck  model  becomes  apparent  at  higher 
current  range.  It’s  necessary  to  examine  their  suitability  and 
accuracy  while  being  applied  in  simulation.  Overall,  LSM  model 
is  found  to  be  the  best  Seebeck  coefficient  model. 

4.  The  internal  temperature  will  be  larger  than  the  hot  side  tem¬ 
perature  beyond  the  reverse  current,  lrev.  The  heat  generation 
inside  will  increase  faster  than  the  cooling  from  the  hot  side. 
Therefore,  lrev  should  be  chosen  as  the  upper  threshold  in  ther¬ 
moelectric  cooler  design  to  prevent  from  high  internal  temper¬ 
ature  which  may  result  in  the  failure  of  test  module. 

5.  Applying  the  CSM  model  is  a  more  conservative  method  to  esti¬ 
mate  the  maximum  temperature  difference  of  TEC.  Actually  the 
temperature  difference  is  larger  in  application. 


6.  The  cooling  performance  can  improve  when  the  Thomson  heat 
maintains  positive.  Joule  heat  dominates  in  the  internal  thermal 
mechanism,  but  more  precise  temperature  control  can  be 
achieved  by  considering  Thomson  heat. 
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